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In this paper we discuss the fluidity of the hot and dense QCD matter created in ultrarelativistic 
heavy ion collisions in comparison with various other fluids, and in particular suggest its possible 
supercriticality. After examining the proper way to compare non-relativistic and relativistic fluids 
from both thermodynamic and hydrodynamic perspectives, we propose a new fluidity measure 
which shows certain universality for a remarkable diversity of critical fluids. We then demonstrate 
that a fluid in its supercritical regime has its fluidity considerably enhanced. This may suggest a 
possible relationship between the seemingly good fluidity of the QCD matter produced in heavy ion 
collisions at center of mass energy of = 200 AGeV and the supercriticality of this matter with 
respect to the Critical-End-Point on the QCD phase diagram. Based on such observation, we predict 
an even better fluidity of the matter to be created in heavy ion collisions at LHC energy and the loss 
of good fluidity at certain lower beam energy. Finally based on our criteria, we analyze the suitabil- 
ity of a hydrodynamic description for the flreball evolution in heavy ion collisions at various energies. 
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I. INTRODUCTION 

The exploration of the QCD phase diagram as well as 
the quantitative characterization of QCD matter is one 
of the most interesting challenges and questions in strong 
interaction physics. Hot and dense QCD matter can be 
created in the laboratory by means of heavy ion collisions, 
and experiments at the CERN Super Proton Synchrotron 
(SPS) and at the Relativistic Heavy Ion Collider (RHIC) 
have, over the years, revealed many intriguing and unex- 
pected properties of this matter. For example the mea- 
surement of an unexpectedly large elliptic anisotropy V2 
[II can be reproduced within the framework of hydro- 
dynamics [2j], at least for low transverse momenta (pt)- 
This observation has led to the conjecture, that the mat- 
ter produced in these collisions is strongly interacting, 
with nearly ideal fluidity j3j. Several microscopic expla- 
nations for this behavior have been suggested |4 1 @ @ , for 
example the "magnetic scenario" which features the co- 
existing electric and magnetic components of the plasma 
with the magnetic one ultimately enforcing the QCD con- 
finement transition. 

Arguments for the nearly ideal fluidity at center of 
mass energies of ^/s — 200 AGeV, quantitatively rep- 
resented by a rather small ratio of shear-viscosity over 
entropy-density, ry/s, have come from different directions. 
Firstly, the V2 data measured at RHIC have provided 
rather stringent constraints on viscous hydrodynamic cal- 
culations: the current status based on this approach is 
that an upper limit rj/s < 6/(47r) may be set [2j. (Some 
caveats on interpretation of V2 data however should be 



kept in mind, such as hidden non-collective contributions 
to V2 or alternative mechanism of producing 1^2 [8]). 
Secondly, recent developments from the AdS/CFT cor- 
respondence [9^ have hinted at a universal lower bound 
v/-^ ^ l/(47r). While in principle violating examples have 
been found (see detailed discussions in e.g. for all 

practical purposes this value serves as a useful bench- 
mark for good fluidity of relativistic fluids. Finally, the 
fluidity of the QCD fluid has been compared [ll|, |T2| to 
a few commonly known substances like Helium-4, water, 
nitrogen etc, using the measure rj/s: it has been con- 
cluded that the QCD fluid has smaller ry/s as compared 
with all those normal substances. 

In this paper we will revisit the concept of fluidity. 
For a sensible discussion of fluidity one needs to clearly 
distinguish the concept oi fluidity from the concept of ap- 
plicability and validity of hydrodynamics. The concept 
of fluidity of a substance can only be meaningful if the 
fluidity is defined exclusively in terms of properties of the 
substance itself. The question of applicability of hydro- 
dynamics, on the other hand, requires the effective mean 
free path^ to be small compared to the size of the sys- 
tem or rather the typical wavelengths of the excitations 
which are supposed to be described within hydrodynam- 
ics. Consequently, in this case we compare a property of 
the system, the mean free path, with an external scale 
characterizing the situation we want to describe within 
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^ To simplify the argument, in tlie rest of the introduction we will 
use the familiar concepts of "mean free path" and particles, not- 
ing that these are limited to a kinetic description of a sufficiently 
dilute system. In section IIIII we will derive a more general ex- 
pression for the relevant length scale which will not require the 
applicability of kinetic theory. In the kinetic limit it reduces to 
the mean free path. 
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hydrodynamics, for example a sound mode of a given 
wavelength. And given a sufficiently large external scale, 
hydrodynamics is always applicable. For example we ob- 
serve sound propagation in both water and air, and yet 
one would be inclined to assign a better fluidity to water 
than to air. Therefore a meaningful definition of flu- 
idity needs to measure the effective mean free path in 
terms of a length scale inherent to the substance itself. 
One obvious choice is the interparticle distance, or more 
generally, the density-density correlation length. This al- 
lows the comparison of substances with vastly different 
inherent length scales, such as water vs. a Quark-Gluon 
Plasma. Obviously, the minimum wavelength (in me- 
ters) for sound propagation even in a weakly interacting 
QGP is still considerably shorter than that of a strongly 
interacting atomic or molecular substance. However, if 
we ask the question "what is the minimal wavelength for 
sound propagation measured in units of the interparticle 
distance?" then a comparison between these substances 
becomes meaningful, and the strongly interacting sub- 
stance will likely exhibit a better fluidity. 

Following these general considerations, in this paper 
we discuss the fluidity of the hot and dense QCD mat- 
ter created at RHIC by comparing it with normal non- 
relativistic fluids and by applying valuable insights from 
those fluids. While such a comparison is not new (see 
e-g- we differ from all previous approaches in a 

few distinct points. First of all, we carefully examine the 
difference between the relativistic and non-relativistic flu- 
ids. This includes their different inertia: the relativistic 
inertia, i.e. the enthalpy, reduces to the mass density in 
the non-relativistic regime. This also includes the choice 
of reference scale, as discussed above: while the temper- 
ature may serve as a reasonable estimator for the inter- 
particle distance for relativistic fluids, it certainly does 
not in the non-relativistic limit. Based on these consid- 
erations, we propose a new measure of fluidity for com- 
paring various fluids. Furthermore we emphasize, for the 
flrst time, the possible relevance of the so-called super- 
critical fluid for the matter produced at RHIC and discuss 
important implications for the expected fluidity for the 
matter produced in heavy ion collisions at different center 
of mass energies. We will also elaborate on potential con- 
sequences for the search of the QCD Critical-End-Point 
via a beam energy scan. 

As just discussed, a closely related, but different, ques- 
tion is the applicability of hydrodynamics in heavy ion 
collisions at various beam energy ^/s. For the descrip- 
tion of the dynamical evolution of system within hydro- 
dynamics to be valid the length scale characterizing the 
variations of the system needs to be large compared to the 
effective mean free path of the (quasi) particles within the 
fluid Thus, good fluidity of the underlying matter 

may not guarantee applicability of hydrodynamics, if the 
typical gradients of the flow field are large, e.g. when the 
system size is very small. Based on an analysis of sound 
wave attenuation, we will provide quantitative criteria 
for the applicability of hydrodynamics, and evaluate the 



situations at SPS, RHIC, and LHC respectively. 

This paper is organized as follows: In Section |lT] we 
will discuss the difference and relation between a non- 
relativistic (NR) fluid and a relativistic(R) fluid. In 
Section IIIII we will then propose a new fluidity mea- 
sure, which is applicable for both relativistic and non- 
relativistic systems. Based on this new measure we will 
compare various fluid systems and demonstrate the im- 
provement of fluidity in a fluid's supercritical regime. In 
Section ITVl we will use Lattice QCD results to construct 
equal-pressure lines on the QCD phase diagram. We then 
discuss the relationship between fluidity and supercriti- 
cality for heavy ion collisions, the evolution of fluidity 
with beam energies, and its implications for the search of 
the QCD Critical-End-Point(CEP). Finally in Section|V] 
the applicability of hydrodynamics in heavy ion collisions 
at various energies will be discussed. 

II. RELATIVISTIC AND NON-RELATIVISTIC 
FLUIDS 

When comparing a relativistic (R) with a non- 
relativistic (NR) fluid, one needs to carefully keep track of 
the mass terms, which customarily are neglected in non- 
relativistic thermodynamics. In non-relativistic thermo- 
dynamics, the basic thermodynamic relation 

ENB^ = TS-pV + ^iNBN (1) 

does not take into account the mass of the particles. Here 
Enb. refers to the kinetic and interaction energy of the 
particles. Similarly, the chemical potential, which repre- 
sents the increase of energy by the addition of one extra 
particle, does not account for the particle's mass. The 
relativistic version of the basic thermodynamic relation, 

ER^TS-pV + ^lRN, (2) 

on the other hand, takes into account the particle masses. 
Its non-relativistic limit can be obtained by simply in- 
cluding the mass terms in both the energy and the chem- 
ical potential 

Er = Enr + m (3) 
^J■R = fJ-NR + rn. (4) 

As we will discuss below (see also [13, [ll]), the thermody- 
namic quantity entering hydrodynamics is the enthalpy 
density, w defined as 

wr^€ + p = Ts + pru (5) 

where e is the energy density, p the pressure, s the 
entropy-density, and n is the particle-density. As we 
will show, the non-relativistic limit of hydrodynamics in- 
volves the non-relativistic limit of the enthalpy, wr, in- 
cluding the mass term 

w = Ts -\- {pnr + m)n mn = p (6) 
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and it is dominated by the mass density. In the uhra- 
relativistic hmit T 3> /^/j, on the other hand, the enthalpy 
density is given by 



Ts. 



(7) 



The kinematic viscosity, which is defined as the ratio of 
the shear viscosity over the enthalpy density 



V 

V = — 

w 



(8) 



usually serves as a measure for dissipation [17|. While 
the widely used ratio of shear- viscosity over entropy den- 
sity, 77/s is indeed related with the kinematic viscosity for 
ultra-relativistic systems, it misses the dominant mass 
term for a non-relativistic fluid. 

In the following, we further demonstrate this point 
both in thermodynamics and in hydrodynamics^. 



Next we examine the non-relativistic limit by taking 
/3mc^ — 00 in function I[y\, EgfTOl (and its derivative 
/' ^ dM/dy): 



f ^ NR: nmc^ + niksT) [in(nA^) - l] 
H NR : mc^ + {kBT)ln{nX^) 



NR : nmc? + —n{kBT) 



(13) 



with A = /mkBTY^'^. Obviously, in the non- 

relativistic regime the energy associated with the rest 
mass dominates the chemical potential, the free energy 
density, the energy density, and the enthalpy density. 
Contrary to that, neither entropy density s — {e — f)/T 
nor the pressure p — /in — /, have an explicit dependence 
of the mass term, as intuitively expected. 



A. Thermodynamics 

We start with the example of a classical (Boltzmann) , 
relativistic free gas at temperature T and fixed (net- 
)particle density n = np — np. Following standard sta- 
tistical mechanics in e.g. [l^, the Helmholtz free energy 
density is given by 



/(r,n)= niksT) 



+ln 




(9) 



with /3 = l/iksT) and A = {2Tr'^)'^/^ {phc). Here we have 
introduced the re-scaled momentum integral, p — l3pc , 



dpp e 



(10) 



Using the usual thermodynamic relations we obtain ex- 
pressions for other quantities, such as the chemical po- 
tential and the energy density: 



dl 
dn 



iksT) In 



■)2+4 



(11) 



= n(kBT)\l 



21 



r 



(12) 



B. Hydrodynamics 

We now turn to the hydrodynamics in the relativis- 
tic and non-relativistic regime. Since the non-relativistic 
limit of relativistic hydrodynamics is discussed in text- 
books [13, [ll], we will be brief here and just remind 
ourselves of the essential points, in particular how the 
relativistic inertia i.e. the enthalpy density w is replaced 
in the non-relativistic limit by the mass density pc^ . We 
start with the non-relativistic Navier-Stokes (N-S) equa- 
tion"^ : 



[dt+v- \/]v = - 



VP ^ V 
P P 



(14) 



with the non-relativistic shear tensor S-'* = djVi + diVj — 
If^jiV ' The corresponding relativistic N-S equation is 



given by^: 



1 



-dop] 



(15) 



with the relativistic shear tensor = c [9^ 
{u ■ d)uf_,u^ + |(u^u^ - gp.^){d ■ u)], 7 = 1/ a/1 - v^/c'^, 
Uf^ = 7(1, w/c), and do — ^dt- In the non-relativistic 
limit one has 7 —J' 1, and w/c^ — p and thus in leading 
order of v/c recovers the non-relativistic Navier-Stokes 
equation |17| . 

To further elaborate the point, let us consider the prop- 
agation and attenuation of a sound wave of frequency w 



^ To better keep track of the mass terms, in the following two 
subsections we make explicit the dependence on the speed of 
light c, the Boltzmann constant kg and the Planck constant h, 
while in the rest of the paper we use natural units with these 
constants taken to be unity. 



^ In this discussion of N-S equation we neglect bulk viscosity and 
assume constant shear viscosity across the fluid. 
It is well-known that for relativistic hydrodynamics the derivative 
expansion to only first order, i.e. the N-S form, has causality 
problem and a consistent treatment requires higher orders in 
derivative expansion, see e.g. 0. 
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and wave vector k = 27r/As (As is the wavelength ) in 
the presence of dissipation, characterized by the shear 
viscosity rj. The dispersion relation for the wave is given 
by 



R fluid 
NR fluid 



(16) 



To take into account possible bulk viscosity ^, one sim- 



ply makes the replacement f?? ^ 
the speed of sound Cs is given by 



[77 + C- Furthermore 




, R fluid 
NR fluid 



(17) 



As one can see, the same correspondence {w, e)^ — > 
{pc^)nr appears again as in the thermodynamics 
Eq.dlSl). 



not the dynamical viscosity r] itself but rather the so- 
called kinematic viscosity rj/ p (see e.g. [l3|)- Indeed by 
looking at actual data for water in the liquid and vapor 
phase at the same pressure one finds that liquid water 
has about one order of magnitude bigger shear viscosity, 
?7, but nonetheless "wins the fiuidity contest" since its 
kinematic viscosity, rj/ p, is about two orders of magni- 
tude smaller than that of the vapor phase. We note that, 
the expression •q/w for a relativistic fiuid is a relativistic 
version of the kinematic viscosity (see also related dis- 
cussions in [3]). To conclude, we emphasize that for a 
non- relativistic fiuid rj/p is different from rj/s (as is evi- 
dent from the thermodynamics discussion) and only rj/p 
serves as a good measure of fluidity (as is evident from 
the hydrodynamics discussion). This also implies that 
using r]/s to compare the fluidity between a relativistic 
fiuid (like the QGP, the AdS/CFT plasma) and a non- 
relativistic fluid (like Helium, water or cold Fermion gas 
[20I ) may not be as informative as one would expect. 



C. Discussion on rj/s 



III. A NEW FLUIDITY MEASURE 



We end this section by a discussion of the ratio of shear- 
viscosity over entropy-density-ratio, rj/s, which has at- 
tracted a considerable interest in various fields of physics. 

We first recall, from the perspective motivated by 
studying the QGP via heavy ion collisions, why rj/s 
is a useful measure for the relativistic fluid, as was 
first pointed out in the seminal paper [l9| . Consider a 
Bjorken-type longitudinal expansion of the quark-gluon 
plasma (with low baryonic density) formed in heavy ion 
collisions: its system size is limited by cr at early time. 
With the presence of shear viscosity rj, there will be dis- 
sipative effect (e.g. entropy generation) characterized by 
the ratio (roughly the Knudsen number) which, by 
thermodynamic relation Eq.Q w ~ Ts, leads to j 
To ensure a controllable dissipative correction to the ideal 
hydrodynamics, one requires firstly rj/s is small enough 
and, secondly, Tt is of the order 1 or bigger. The former 
shall be a property of the underlying QGP while the lat- 
ter means the viscous effect is most severe at early time 
and the hydrodynamic evolution may not be a good ap- 
proximation at too early time. (Note in this discussion 
we adopted natural relativistic units with c = 1). One 
can draw two conclusions from this discussion: (a) rj/s 
can serve as a good measure of fluidity for a relativistic 
fluid and the smaller it is the better the fiuidity; (b) the 
ability of 77/s to serve such a role is actually inherited 
from r]/w. 

The discussion above, however, leads to the observa- 
tion that 77/s for a non-relativistic fluid does NOT nec- 
essarily provide a good measure for its fluidity, because 
the role of the relativistic rj/w corresponds to the non- 
relativistic 77/p. This is certainly not a new lesson: it has 
been known from Navier and Stokes's time that what 
really matters for usual non-relativistic fiuids' fiuidity is 



In this section we propose a new fiuidity measure that 
is suitable for comparison between relativistic and non- 
relativistic fiuids. In the following we will analyze the 
propagation of sound modes in the presence of dissipa- 
tion (viscosity) and ask ourselves under what condition 
the dissipative (viscous) terms in the equations prevent 
sound from propagating. The advantage of this strategy 
is that we stay entirely within the framework of viscous 
hydrodynamics and do not need to make additional as- 
sumptions, such as the applicability of kinetic theory. As 
a consequence our fluidity measure will be expressed in 
terms of well defined quantities such as the shear viscos- 
ity, the speed of sound, etc. 

We start with the sound dispersion relation in Eq. ([T5| . 
By requiring the imaginary part of the frequency, Imuj, 
to be small in magnitude as compared with its real part 
TZeuj, we obtain: 



TZeuj 



« 1 
A. 



27r 


47r 

T 


r 




1 (to/c2 
1 — 


)Cs 







R fiuid 



(18) 
(19) 



The above equation implies that if a sound wave has its 
wavelength Xs comparable (or even smaller) than Lrj, it 
will be quickly damped on (or shorter than) a time scale 
of its period and spatially on a length scale about (or 
shorter than) its wavelength, which essentially means it 
can not propagate away in the medium. Therefore, the 
physical meaning of the length L,j introduced above is 
to provide a measure for the minimal wavelength of a 
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sound wave to propagate in such a viscous fiuid^. A 
more quantitative criteria will be discussed in Section IVl 
and given in Eq.([27|). 

Furthermore the length has the meaning of an effec- 
tive mean-free-path (MFP) in terms of microscopic fluid 
particle motion. This becomes transparent in a weakly 
coupled gas: taking the non-relativistic gas as an exam- 
ple, the shear viscosity according to kinetic transport is 

77 pvtImfp (20) 
while the speed of sound is 

with vt the thermal velocity. These lead to the combi- 
nation 

Ln — — ^ ^ Imfp (22) 

PCs 

Unlike the mean-free-path which is conceptually intuitive 
but practically not easily computable or measurable (e.g. 
for fluids), the length is well-defined by macroscopic 
properties of the fluid and thus of practical use. 

We emphasize that involving the speed of sound in the 
definition of length scale is essential. This can be seen 
already from a dimensional argument: r]/ p or ri/{wlc^) 
has the dimension of [Length]^ [Time]^^. To turn this 
into a length scale one needs to divide it by a quantity 
of the dimension [Length] [Time]~^ i.e. a velocity. The 
natural choice here is the speed of sound as a charac- 
teristic of the macroscopic matter. The speed of light, 
c, on the other hand is not suitable here since it is nei- 
ther a property of any specific substance nor should it be 
of any relevance to non-relativistic fluids, such as water. 
Since the speed of sound changes considerably close to a 
phase-transition or a rapid crossover, its inclusion in the 
fluidity measure J- and in the effective mean free path 
is essential. 

Next we need to introduce another meaningful length 
scale to make a dimensionless ratio: this becomes a ne- 
cessity when comparing fluids at vastly different scales. 
In case of the well known dimensionless ratios like the 
Knudsen and Reynolds numbers, an external length scale 
characteristic of the fluid motion is introduced, like the 



^ We note that other types of dissipative processes like ther- 
mal conduction may also be present and thus introduce differ- 
ent length scales. For example when the sound wave period 
Ts ~ l/(csfe) becomes larger than the thermal relaxation time 
scale TT ~ l/{DTk^) (with Dt the thermal diffusivity) at wave- 
length smaller than Dx /cs , then heat transport becomes rather 
efficient and the sound propagation becomes isothermal (see ex- 
amples in e.g. \2i ). Nevertheless additional sources for dissipa- 
tion do not change the fact that the "good" sound modes shall 
have their wavelengths (at least) larger that the Lrj set by shear 
viscosity only. 



diameter of a pipe or the size of a moving object inside 
the fluid, etc. In the context of relativistic fluids created 
in heavy ion collisions, there have been scaling studies 
of collective flow for lower energy collisions based on the 
Reynolds number in [l3| and more recently for higher en- 
ergy collisions based on the Knudsen number in |I4| | (see 
also related discussions in [l^ ) . The external length scale 
in those numbers, however, is not an intrinsic property of 
the fluid that we would like to invoke for comparing flu- 
ids across vastly different scales. Instead, a de-correlation 
length of certain density-density spatial correlator gives 
a natural scale of short-range order in the system. In 
most cases this de-correlation length is simply set by the 
inter-(quasi-)particle distance. For a non-relativistic fluid 
"particles" and their number density n are well-defined, 
and thus the inter-particle distance is also a well-defined 
length scale: 

Ln = (23) 

713 

For relativistic fluid it is less straightforward. For ex- 
ample a QGP with ps = and thus ub = can still 
have substantial numbers of quarks and gluons. For 
such a relativistic fluid a simple estimate can be made 
via the entropy density, i.e. n ~ 4§~- Another way 
is to calculate the de-correlation lengtfi of e.g. the cor- 
relator < Tno(x, 0)Tnn(0, 0) > as has been done in re- 
cent lattice work [22| which finds a short range order 
about 0.6/T for QGP in I — 2Tc, in reasonable agree- 
ment with estimate from n -jf—- There could still 

4fcB 

be academic examples where the entropy density can 
be infinite while the short-range order does not van- 
ish: AdS/CFT offers such an example in which the en- 
tropy density goes as s oc N^T^ — > 00 (in the large 
Nc limit) implying — > while spatial correlators 
like < T^y{x^ 0)T^i/(0, 0) > give a non-zero de-correlation 
length ^ l/T. We will return to these issues in the Sub- 
section nil CI 

Finally by taking a ratio of the two length scales, we 
arrive at a fluidity measure 

. (24) 

Below we will flrst show that the measure works well for 
non-relativistic fluids and bears certain universality for 
critical fluids. We will then show that the so-called su- 
percritical fluids have even better fluidity. At the end we 
present a comparison of various interesting fluids. Fol- 
lowing [ll|, , we have extensively used the measured 
data for various fluids from the NIST WebBook [H]. 

A. Critical Fluids 

We first examine various fluids at fixed critical pressure 
P = Pc- In FiglUleft), the proposed fluidity measure 
F — Ljj/Ln is plotted for fifteen different substances at 
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FIG. 1: (Color online) (left panel) Fluidity measure T — Ln/L„ versus T/Tc for fifteen different substances at fixed critical 
pressure P = Pc, see text for more details. The sharp peaks centered at Tc are due to the actual divergence of the shear 
viscosity at the critical point, (right panel) The minimum of the fluidity measure for various substances (as obtained from 
respective curves in the left panel) versus their molar masses, with the stars (from left to right) for H2, ^-ffe, H2O, D2O, Ne, 
N2, O2, Ar, CO2, Kr, Xe, and the diamonds(from left to right) for C4H10, Cg-ffis, C12H2G, C^Fs- 



their respective critical pressure including: Hydrogen 
(i?2), Helium-4 {^He), Water {H2O), Deuterium oxide 
(D2O), Neon (iVe), Nitrogen (iV2), Oxygen (O2), Ar- 
gon {Ar), Carbon Dioxide (CO2), Krypton [Kr), Xenon 
(Xe), Isobutane [C^Hiq), Octane (Cg^fig), Dodecane 
(Ci2-ff26), Octafluorocyclobutane {C^F^). These sub- 
stances cover a wide range of molar mass, chemical struc- 
ture and complexity, with their respective critical tem- 
perature Tc and pressure Pc differing by orders of magni- 
tude. Despite such huge differences, their fluidity curves 
resemble each other not only in shape but even quan- 
titatively. In particular in their good liquid regime — 
roughly the "valley" region at ^ 0.7 — ITc — they 
all show amazingly similar fluidity. To further expose 
the similarity, in the right panel of FiglT] we show the 
value of the fluidity measure J- at its minimum versus 
the substances' molecular molar masses. From this plot, 
roughly two bands can be identified: the green stars 
spread in a narrow band of G (0.3,0.45) while span- 
ning two-orders-of magnitude in molar mass which in- 
clude all 11 non-organic substances; the red diamonds 
with roughly twice bigger F , include 4 organic substances 
with much more complicated molecular structures (e.g. 
chains) which, not surprisingly, lead to more dissipation. 
Even so, the splitting in fluidity between the two bands is 
merely a 0(1) factor rather than any order-of-magnitude 
difference. Nevertheless as one can imagine, with increas- 
ing chemical complexity and molecular mass, the fluidity 
of more complex systems like e.g. engine oil may deviate 
significantly from what are shown in the figure. 



To conclude the study of critical fluids, there appears 
to be certain universality of the newly proposed fluidity 
measure F, indicating that "a good fluid is a good fluid" 
despite many other details regarding the microscopic de- 
grees of freedom. 



B. Supercritical Fluids 

While much attention has been paid to critical fluids, 
the fluidity of a fluid with significantly larger pressure 
P » Pc has been little discussed in the context of heavy 
ion collisions and QCD matter: so let us next explore 
this region. In the literature, this region is often referred 
to as the "supercritical fluid" region 21], deflned on a 
typical substance's T — P phase diagram (with critical 
point Tc, Pc) as: 



supercritical : T > Tc &z P > Pc 



(25) 



In particular we want to explore the deeply supercritical 
region with P » Pc- Taking water as an example, in 
FigHJleft) we plot the fluidity measure J" versus T/Tc for 
fifteen different fixed pressure values ranging from 5 MPa 
ah the way to 1000 MPa (note for water Pc = 22 MPa). 
The dashed (green) curve is for fixed critical pressure 
P = Pc = 22 MPa, while the four solid (red) curves 
above it are for P < Pc, and the ten solid (blue) curves 
below it are for P > Pc- As can be seen from the plot, 
the curves change shape gradually and the fluidity be- 
comes better and better with increasing P. The "valley" 
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FIG. 2: (Color online) The dimensionless fluidity measure T (left panel) and the length scale, L,, (in units of A) (right 
panel) versus T/Tc for water at flfteen different flxed pressure values. In both panels, the dashed (green) curve is for 
fixed critical pressure P — Pc — 22MPa, while the four solid (red) curves above it are for P < Pc, and the ten 
solid (blue) curves below it are for P > Pc, with the respective pressure values for each curve (from top down) being 
P = 5, 10, 15, 20, 22, 25, 30, 40, 50, 100, 150, 200, 250, 500, 1000 MPa. 
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FIG. 3: (Color online) The shear-viscosity-entropy-density- 
ratio rj/s (in unit ft/(47rfc_B)) versus T/Tc for water at flfteen 
different fixed pressure values. The dashed (green) curve 
is for flxed critical pressure P = Pc = 22 MPa, while the 
four solid (red) curves above it are for P < Pc, and the ten 
solid (blue) curves below it are for P > Pc, with the respec- 
tive pressure values for each curve (from top down) being P — 
5, 10, 15, 20, 22, 25, 30, 40, 50, 100, 150, 200, 250, 500, 1000 MPa. 



where J- remains small and relatively flat becomes much 
wider and eventually flattens out substantially above Tc 
for P >> Pc. To quantify the remarkable fluidity of 
the supercritical fluid, we note that the minimum on the 
P = Pc curve has a value for the fluidity Frain{Pc) ~ 0.33 
while the minimum for the P = 1000 MPa ~ 45Pc curve 
is at J^mi„(45Pc) w 0.11, getting smaller by a factor of 
3 and remaining close to the minimum within a rather 
broad temperature region! 

It is also interesting to examine whether the length 
scale Lr^ itself shows similar trends, 3jS L'q is the essential 
scale for discussing the applicability of hydrodynamics 
where it is to be compared with the external scale char- 
acterizing the variation of flow fleld. In the right panel of 
Fig 12] we plot L,, itself for the same conditions and flnd 
that, similar to the fluidity measure T , the scale L,, also 
becomes considerably smaller in supercritical water. 

It should be mentioned that the same observation is 
also true for other fluids that we examined, like IIelium-4, 
Nitrogen, etc. Furthermore such behavior is not specific 
to our fluidity measure. In Fig 13] we plot the widely used 
ratio of shear viscosity over entropy-density, 77/5, which 
show the same qualitative behavior. 

The main lessons, as we emphasize again, are (a) for 
a given substance, the best fluidity is not necessarily 
achieved close to the critical point and (b) when going 
deeper into the supercritical regime its fluidity becomes 
much better. 
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FIG. 4: (Color online) Comparison of fluidity measure for 
various fluids (see text for more details). The curves with 
question marks indicate current estimates of the respective 
fluidity with possible uncertainty, while the curves for Helium 
at Pc and for water at Pc, llPc,45Pc are from actual data. 



C. Comparison of Various Fluids 

Finally we attempt to compare various fluids of great 
current interests in terms of the nevi^ fluidity measure T 
we have proposed and studied above. The results are 
shown in FiglH Below we give the details of how the 
curves for QCD, Cold Fermi Atom, and AdS/CFT are 
obtained. 

For the QCD system, we have used a parametrization 
of the viscosity 77 by Hirano and Gyulassy in [^J]: one 
uses 'q/T^ w T/T^ for hadronic gas (H.G.) below T^, and 
7y/T3 « (T/Tc)3[l + yV(T)ln(r/rc)]2 for the quark-gluon 
plasma(sQGP) in the temperature interval T/Tc G [1 — 3] 
with 

[W(T)/(47r)] = (9/32) . [SQiT^KsB\n{An / g\T))]-^ 
interpolating to pQCD results for T >> Tc, where the 
parameters are given as /3o = 10, Ksb = 12 and the 
running coupling 

[g^iT)]-^ = (9/8^2)in(2^r/A) + 

(4/97r2)ln(21n(27rr/A)) (with A w 190MeV) ( see 
[2J][ll| for more details). The enthalpy density 
w = e + p and speed of sound c, are taken from recent 
lattice results by Karsch et al [25| for 2+1 flavor QCD 
with ~ 22{)MeV . As we mentioned before, L„ is 
estimated byl/(s/4fcs)^^^ with the entropy density also 
taken from j25| . 

For the strongly coupled AdS/CFT system, the shear 
viscosity is well known to be 77/s = l/(47r) [9]. As we 
also pointed out before, there is a short-range order at 
the length scale L„ ~ 1/T however the pre- factor is not 



accurately determined. We simply use L„ = 1/T as an 
estimate. This gives the fluidity T = V3/(47r) w 0.138. 

For the Cold Fermi Atom gas, its shear viscosity has 
been measured near its Feshbach resonance by Thomas 
et al in (26| for a certain range of system energy by an- 
alyzing the damping of collective modes in the atomic 
cloud (see also related work in [l^]). As a benchmark, 
we take the lowest viscosity found from the measurement 
(see Fig. 4 in [2^) which is « 0.214 ftn. The speed of 
sound has also been measured by the same group in psj , 
from which we take the value Cs ~ 0.3632 ti/^ near the Fes- 
hbach resonance. Since the Fermi velocity vp = hkp/m, 
the mass density p = mn, and the Fermi momentum 
kp/n^^^ = {diii^Y/^, we obtain for the fluidity measure 
J- w 0.214/(0.3632 • (37r2)i/3) ^ o.l91. 

The question marks in Fig|4] indicate that the current 
estimates presented above may carry sizable uncertain- 
ties, and we expect the knowledge on these systems will 
become more accurate with time. 

We finally come to the comparison in Fig|4l As one 
can see, the critical fluids (water and Helium-4 at Pc) 
have a somewhat worse fluidity than the QCD, AdS/CFT 
and Cold Fermi Atom systems. Supercritical water at 
P — 11 Pc, on the other hand, already has a fluidity 
comparable to the QGP while at P = 45Pc the fluidity 
of supercritical water appears to be better than that of 
the recently discussed "nearly perfect fluids" . 



IV. A SUPERCRITICAL QCD FLUID AT RHIC? 

The fluidity study in the previous section, in particular 
on the supercritical fluid, naturally leads to the follow- 
ing interesting question: Are we observing a supercritical 
QCD fluid at RHIC? 

As is well known, what is usual referred to as Tc ~ 
170 ^ 190MeV in QCD is not a true second-order phase 
transition temperature but rather the temperature for a 
rapid crossover at ^ = [l^. There is, however, a (hypo- 
thetical) Critical-End-Point(CEP) at (Tcep, t^CEp) on 
the QCD phase diagram marking the end of a plausible 
first-order phase transition at low T but high fi (see e.g. 
[29I I and references therein). Accurate determination of 
the CEP from lattice QCD [30l - l33 | and unambiguous ob- 
servation of the CEP from heavy ion collisions [l2l.f35ll36j 
are among the most interesting and exciting goals of QCD 
research. While currently the position of the CEP is not 
well constrained, it is expected to have lower temperature 
than that of the QCD crossover transition at vanishing 
baryon density, Tcep < Tc- Being very aware on the 
present uncertainty of the actual location of the QCD 
critical end point, let us, solely for demonstration pur- 
poses and for the sake of the argument, assume that its 
location is close to the estimate of ref. (s^l , which sug- 
gests that Tcep ~ OMTc^icep « I.STcep- This loca- 
tion is indicated in FiglS] as the filled black circle along 
with a question mark. 

Next if we adopt the definition of supercriticality as in 
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FIG. 5: (Color online) Schematic isobaric contours (i.e. with 
constant pressure along each line) on the QCD T — /i phase 
diagram with the filled black circle indicating a possible posi- 
tion of the hypothetical Critical-End-Point (CEP); the dashed 
blue horizontal short line stretching from the CEP to the right 
is included to indicate the T = Tcep boundary (see text for 
more details). The filled red, blue and green circles indicate 
estimates of the initial (T, /x) reachable at SPS, RHIC and 
LHC, respectively. 



Eq. ([25|) for QCD, we need to know the constant-pressure 
lines on the QCD phase diagram typically plotted in 
terms of T — /i. To schematically construct these lines, we 
make use of the Taylor expansion of the pressure P{T, fi) 
with respect to /i: 



24 



xf(T) 



(26) 



with X2n = di u/TY ^ " I ^—0 baryonic susceptibilities 
[3^ [33I [37I Issl . l40j . In order to construct the equal- 
pressure lines shown in FigjS] we have used the lattice 
gD resuhs from Cheng et al for Po{T) 25] and xf , xf 

The blue curve in FiglS] going right through the indi- 
cated CEP point indicates the constant-pressure line with 
fixed critical pressure P — P{CEP) which, together with 
the T — Tcep line (dashed blue horizontal one stretch- 
ing from the CEP to the right) form the boundary above 
which there is the QCD supercritical region. We empha- 
size that above the boundary the pressure P ^ T'^ in- 
creases very rapidly. As a result the QGP quickly enters 
deeply into the supercritical regime where P ^ P{CEP). 
To give an idea: the pressures of the lines in units of the 



critical pressure P{CEP) (value on the blue curve) are 
(from bottom to top) 0.05, 0.08, 0.13, 0.23, 0.27, 0.35, 
0.63, 1, 2.7, 10, 18, respectively. Of course with improved 
lattice results for pressure, CEP and (higher order) sus- 
ceptibilities a more accurate equal-pressure map for QCD 
matter will become available. 



We now discuss a few implications for heavy ion colli- 
sions experiments. As a precaution, the following discus- 
sions are by no means intended to provide quantitative 
statements but rather qualitative yet interesting ideas. 

(I) At current RHIC energy = 200 AGeV), FiglS] in- 
dicates that most of the QGP phase of the created mat- 
ter is likely in the supercritical region. This may very 
well be the reason why such good fluidity has been ob- 
served. If a relationship between good fluidity and su- 
percriticality is indeed true for the QCD matter just as 
it is for water, then there will be a sensitive dependence 
of the fluidity on the actual position of the CEP. For 
example, the matter created at RHIC might have most 
of its evolution being in the supercritical region in case 
Tcep and Pcep £^re considerably lower than indicated 
in the FiglS] On the other hand, the matter created at 
SPS [y/s ~ 20 GeV) seems to have its initial pressure 
and temperature already rather close to the critical one 
which would result in poor fluidity. In principle these re- 
gions (supercritical, near-CEP, below-CEP) of very dif- 
ferent fluidity can be accessed experimentally by tuning 
the center of mass energy ^/s which changes the initial 
densities (see e^. [4l|) as well as the freeze-out points 
(see e.g. [12] [13) in such collisions. To which extend 
this difference in fluidity transforms into measurable ef- 
fects, such as non-hydrodynamic behavior and thus less 
elliptic flow, is a non-trivial question. As we will discuss 
in the next section, besides the fluidity the actual size 
of the system is essential for the applicability of hydro- 
dynamics. Our estimates (see FiglHI) show that even for 
the very good fluidity assumed for a strongly interacting 
QGP, denoted as "sQGP" in FigH] a hydrodynamic de- 
scription for the flreball expansion is a best marginal at 
SPS energies. The reason is that, according to our esti- 
mate, the system starts close to the phase transition and, 
thus, reaches the hadronic phase while still comparatively 
small in size. Consequently, it may be difficult to tell if 
the system enters a regime of reduced fluidity from flow 
studies alone. Of course if our estimates are wrong and 
one finds evidence that the initial entropy density reached 
is considerably larger, a change in the fluidity may show 
up in the systematics of elliptic flow measurements, such 
as its dependence on beam energy and system size. 

(II) Based on the possibility of a mostly supercritical 
QCD matter to be created in experiments at the Large 
Hadron Collider (LHC) (see FiglSj), one is led to pre- 
dict an even better fluidity to he observed at LHC than 
at RHIC. And since for the same fiuidity the conditions 
for hydrodynamics are more favorable for LHC energies 
than for RHIC, we would expect nearly ideal hydrody- 
namic evolution at LHC. 
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V. APPLICABILITY OF HYDRODYNAMICS IN 
HEAVY ION COLLISIONS 

In this section we turn to a discussion of applicability 
of hydrodynamics in heavy ion collisions, particularly for 
SPS, RHIC and LHC. As already discussed in the intro- 
duction, the fluidity of the underlying matter is not the 
only determining factor: the "best" (in terms of fluidity) 
fluid may cease to flow if put into a sufficiently narrow 
pipe. In other words, the applicability of hydrodynamics 
depends equally on both the internal properties of the 
substance, i.e. its fluidity or L^, and the external set- 
tings, i.e the typical length scale Ls of the variation of the 
flow field in the system under consideration ^ : one useful 
criteria, equivalent to the famous Knudsen number, is a 
ratio of the two Ljj/Ls. 

In order to derive a quantitative criterion, we again 
use the example of sound dispersion discussed before. By 
requiring that a propagating sound wave has to complete 
at least one full period before its amplitude damps by 
a factor due to the imaginary part, we arrive at a 

Q 2 

minimal wave length Xmin = ^~Lri- For a sound wave 
the length scale characterizing the flow fleld variations 
simply the wavelength, Ls ~ A. Thus we arrive at the 
following criterion for the applicability of hydrodynamics: 



The length L^{T) for QCD (at = 0), shown in 
FigHJleft), is determined as described section ITlI CI It 
relies on a parametrization for rj{T) from the Hirano- 
Gyulassy [2J1 , and we note that the actual shear viscosity 
of QCD matter, once it is known, may be quite different. 
The temperature dependence of shown in the plot 
features a steep drop from the hadronic phase to about 
T ~ l.lTc, followed by a rather constant value for tem- 
peratures above, T > I.IT^. One may expect becomes 
large again at very high temperatures T » Tc where a 
weakly coupled QGP describable by pQCD takes over. In 
principle also depends on the baryon-number chemi- 
cal potential, jis, and this dependence shall be taken into 
account when discussing collisions at low energy (such as 
SPS) where could be sizable. However currently very 
little is known about such dependence, and we here use 
the zero ^ result also for SPS. 

We next give a rough estimate for the typical length 
scale Ls of flow field variance for heavy ion collisions at 
different ^/s ^. At early times, right after the collision. 



the system size is limited mainly by its longitudinal ex- 
tent and the scale governing the variation in the fluid can 
be estimated to be about 2t. At late time when t be- 
comes larger than the nuclear radius ~ Ra the limiting 
scale is set by the transverse size of the fireball which is 
slightly larger than Ra and slowly grows. For the purpose 
of our rough estimate, we adopt a simple prescription: 
Ls ^2TfoTT <R and Ls « 2i? for r > i? with R = 8 fm 
for AuAu and PbPb collisions. For the collision dynam- 
ics which determines the temperature evolution s(r) (or 
inversely t(s)), we simply use the Bjorken flow relation 
s(r) = so{tq) ■ tq/t, which up to r '-^ i? is a reasonable 
approximation. For the initial condition we follow esti- 
mates used in hydrodynamic modelling (see e.g. [4l[) and 
assume equilibration at tq = 1 fm, with flreball center en- 
tropy densities (in central collisions) to be sq = 24 fm~'^ 
at SPS, So = 70fm-3 at RHIC, and sq = 154 fm'^ at 
LHC. 

Finally we turn to the discussion of the applicability 
of hydrodynamics given the criteria of Eq l27l The ratio 
Ln/Ls in units of 3/(87r^) is plotted in FiglHK right) for 
SPS (red diamonds), RHIC (blue stars), and LHC (green 
boxes). The thick dashed horizontal line indicates the 
equality of the condition expressed by Eql^Tl If the sys- 
tem is considerably below this line a hydrodynamic de- 
scription should be a reasonable approximation. If it 
flnds itself considerably above corrections to hydrody- 
namics from higher orders in derivatives or an even non- 
hydrodynamic description are called for. To quantify the 
previous statement, if the system finds itself at a value 
of (L^/L^)/(3/87r2) = 2, such as the points for SPS 
energies on FigjSKleft), then the amplitude of a sound 
mode after propagating the distance of one wavelength 
will be reduced by a factor of ~ 7, and for a value of 
(i^/L,)/(3/87r2) = 5 the reduction would be ~ 150! 
Given our estimate and the assumptions which it is based 
on, a hydrodynamic description is not a good approxi- 
mation for SPS, whereas for RHIC and even more so for 
LHC energies, it seems to be more or less justified*. Such 
differences from SPS to RHIC and LHC lie in the differ- 
ent time evolution due to the different initial densities: 
when cooling down into the region where starts ris- 
ing abruptly near and below Tc, the SPS fireball still has 
a rather small size (longitudinally) while the RHIC and 
LHC fireballs already becomes large with a size about 2i? 
and hence have their hydrodynamic evolution extended 
much longer. While our estimates contain uncertainties 
and should by no means considered to be precise, we have 



We re-emphasizc that the length scale L,, is a well defined quan- 
tity depending only on macroscopic variables that are measurable 
and/or calculable. Furthermore unlike the mean-free-path which 
depends on quasi-particle picture (see e.g. discussions in [Tsll ). 
the length scale L,, remains a useful and meaningful scale even 
for strongly coupled systems. 

In this discussion we focus only on central collisions which make 



the estimation easier. Generally speaking when going from cen- 
tral to peripheral collisions one expects hydrodynamics to be less 
and less applicable. 

We note that even for RHIC and LHC the first fm or so of the 
evolution is above the "criteria line" due to the rather small sys- 
tem size. Whether this has a profound effect (such as in entropy 
generation) or not is unclear at this time (see a very interesting 
discussion in [44|'). 
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FIG. 6: (Color online) (left panel) The length scale Lr, — ^^^^i-^^ (in unit fm) versus T/Tc for QCD at = 0, with r\ from 
parametrization in [23| and w^Cs from lattice data [2^1 . The dashed horizontal lines at 0.1 fm and 0.5 fm are included to guide 
the eyes, (right panel) The criteria (|27|) for applicability of hydrodynamics versus T/Tc at SPS(red diamonds), RHIC(blue 
stars), and LHC (green boxes), with the thick dashed horizontal line indicating the borderline below which hydrodynamics may 
be a good description according the criteria (see text for more details). A line for CuCu(orange triangles) collisions at RHIC 
energy (200AGeV) is also presented. 



provided a semi-quantitative picture for evaluating the 
applicability of hydrodynamics with varying -^5, which 
implies at LHC the fireball expansion shall be even better 
described by hydrodynamics as compared to RHIC. 

Besides the beam energy a/s, a change of the system 
size should also affect the conditions for the applicability 
of hydrodynamics. For example, RHIC has done exper- 
iments with both AuAu collisions and CuCu collisions, 
both at full energy, and it would be interesting to examine 
and compare the applicability of hydrodynamics for these 
two different systems. In Fig 15] (right) we have included 
a calculation of the same applicability measure for CuCu 
(orange triangle) collisions at 200AGeV, with a initial 
center entropy density sq = 45 fm~'^ and a transverse 
size parameter i? = 5 fm. The curve indicates that hy- 
drodynamics is marginally applicable for such collisions, 
as both the initial density and the transverse size of the 
matter formed in CuCu collisions are smaller that those 
in AuAu collisions at the same energy. 



VI. SUMMARY 

In summary, we have discussed the fluidity of the hot 
and dense QCD matter as produced in ultrarelativistic 
heavy ion collisions in comparison with various other, 
well known fluids. In particular we have suggested its 
possible supercriticality based on insights gained from 
studying more conventional fluids like water. We have 
discussed several aspects relevant for a proper compari- 
son of non-relativistic and relativistic fluids, both from 



thermodynamics and hydrodynamics perspectives. A 
new fluidity measure T = L^^/ Ln is then proposed, which 
shows certain universality in the good fluid regime for a 
remarkable diverse set of critical fluids. We have further 
demonstrated that the fluidity is enhanced in the super- 
critical fluid regime on a fluid's phase diagram. These 
studies inspired us to conjecture that the seemingly good 
fluidity of the QCD matter at RHIC may actually be re- 
lated to its supercriticality with respect to the Critical- 
End-Point on the QCD phase diagram. This observation, 
if true, has far-reaching consequences for heavy ion col- 
lisions experiments: (a) the loss of such good fluidity at 
certain lower beam energy which is sensitively related to 
the position of the long sought CEP; (b) an even better 
fluidity may be expected at higher beam energy, which 
will soon be tested by the LHC heavy ion program. Fi- 
nally we have analyzed the applicability of hydrodynam- 
ics for the fireball evolution in heavy ion collisions at 
various energies. Our analysis was based on a model 
parametrization for the shear- viscosity and a rough esti- 
mate of the relevant length scales governing the fireball 
expansion. Given these assumptions we find that hydro- 
dynamics should be applicable for collisions at RHIC and 
LHC energies but not for SPS energies. 
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